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s ■ 

Determination of the value of the leptonic CP-violating phase 8 and the neutrino mass 
hierarchy sgn<5m3! through long baseline neutrino oscillation experiments is systematically 
analyzed. We note that the two oscillation spectra are difficult to discriminate and lead to the 
degeneracy when they are peaked at the same energy and have the same peak probability. 
q ' The condition of peak-matching is therefore introduced as a criterion for the presence of 

degeneracy. The matching of peaks is visualized as an intersection of trajectories traced by 
the peak of an oscillation spectrum while the value of 5 is varied from to 2n. We numerically 
' calculate a pair of trajectories for a pair of hierarchies and examine the degeneracy, especially 

\ that concerning the hierarchy. We formulate the trajectory in terms of analytic expressions 

and evaluate the critical length, which is shown to be proportional to 1/ sin #13. In view of 
■ our analysis, we discuss future prospects to solve the hierarchy degeneracy with regard to the 

' following four approaches: elongating the baseline length sufficiently, using both neutrinos 

and anti-neutrinos, combining experiments with different baseline lengths, and observing two 
. or more oscillation peaks. 

o ■ 

m : , . 

<^> §1. Introduction 

t> ■ 

Experiments suggest the existence of oscillation among the different flavors of 
I neutrinos, providing rich information on the flavor structure of the lepton sector. 1 ) -4 ) 

i" 1 ■ Despite this progress, however, present knowledge of neutrinos is still incomplete. In 

Ph. the notation of Ref. 5), the values of the mass parameters and the mixing parameters 

yet to be determined include one of the mixing angles #13, the sign of <5m 2 i, and the 
CP-violating phase 5. The value of #13 is confined only by an upper bound as 
sin 26*13 < 0.19, 5 ) while the sign of Sm^ and the value of S are completely unknown 
to date. CP violation manifests itself only in the flavor-changing appearance channel 
such as Vn — > u e , which in turn is suppressed by a small factor of sin 2 26*13. For 
this reason, we expect a two-staged program in pursuing the unknown properties of 
neutrinos. The first stage is the search for #13. It is anticipated that its upper bound 
will be improved down to sin 2 2#i3 < O(10~ 2 ) by planned experiments using nuclear 
reactors 6 ^ and accelerators. 7 ) The second stage is the search for the sign of Sm^ 
and for CP violation. The next generations of long baseline neutrino oscillation 
experiments offer promising opportunities for such searches, 7 ) and other probes of 
the sign of <5m 2 i are also discussed in the literature. 8 ) 

Our focus in this paper is on the second stage based on the optimistic expectation 
that the next generation of reactor neutrino experiments will find evidence of non- 
vanishing #13 and that its value will turn out to be large enough for the CP-violation 
search. We consider the search for the CP-violating phase and the sign of Sm^ 
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through long baseline neutrino oscillation experiments using a conventional beam of 
muon neutrinos. A single long-baseline experiment, however, does not necessarily 
lead to the determination of a unique set of values of oscillation parameters. One 
obstacle is the ambiguity of parameter values due to the experimental errors. 9 )' 10 ) 
Another is the parameter degeneracy, which arises when the experimental results can 
be attributed to two or more sets of parameter values. Three types of the degeneracy 
are customarily referred to: the intrinsic degeneracy, on (#13, (Jcp); 11 ^ the hierarchy 
degeneracy, on Sm^ ^ 0; 12 ) and the octant degeneracy, on #23 ^ 7r/4. 13 ) Although 
the parameter degeneracy should be avoided, n )~ 16 ) its presence is difficult to predict 
comprehensively due to the complicated dependence of experimental results on many 
parameters and on various setups of experiments. The plot introduced in Ref. 12) 
gives an overview of the presence of degeneracy and has been found versatile for its 
analysis. It presents the two oscillation probabilities or event rates of two channels 
in a two-dimensional space, enabling to show separately the CP-violating effect and 
the matter effect. 

We introduce another aspect of the analysis of the energy spectrum of the ap- 
pearance probability. Giving an intuitive illustration of the determination of the 
parameters from this spectrum, we offer a view of ours on the emergence of de- 
generacy and its resolution. The pivot of our study is the peak of the oscillation 
probability, particularly its position, or the energy and probability at the peak. We 
note that two oscillation spectra whose peak positions coincide are difficult to distin- 
guish and likely to cause degeneracy. We trace the peak position varying the values 
of 5 and sgn 5m\i to show how their search is entangled and the hierarchy degeneracy 
results. We vary the baseline length as well and describe from our point of view how 
the degeneracy disappears when the baseline becomes sufficiently long. We organize 
these analyses by deriving analytic expressions of the oscillation probability at the 
peaks and offer an overview concerning the presence of degeneracy and possible ways 
to avoid it. 

The outline of this paper is as follows. In £j2j we introduce the peak- matching 
condition as a criterion for the presence of degeneracy and facilitate its visualization 
by drawing closed trajectories of the oscillation peak. In we develop an under- 
standing of the presence and absence of the degeneracy of parameters through the 
loops of the peak. In SJH we derive an analytic expression of the — ► u e appearance 
probability at the oscillation peaks to elucidate its dependence on the mass param- 
eters and the mixing parameters, and show how the loops are distorted due to the 
change of the baseline length. In £j5l we apply the peak loops to the evaluation of 
four methods to uniquely determine the value of the CP-violating phase and the 
mass hierarchy. Section [6] presents the conclusion and discussion. 

§2. Peak-matching condition and the peak loop 

We assume that there are three neutrino generations and adopt the definitions 
given in Ref. 5) of the quadratic mass differences Smfj C {1,2,3}), the mix- 

ing angles 9ij, and the CP-violating phase 5. All experimental results obtained to 
this time can be accounted for by neutrino oscillation by taking <5m|i — (8.OI03) x 
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10" 5 eV 2 , {Sm^l ~ (1.9-3.0) x 10~ 3 eV 2 , sin 2 26> i2 = 0.86+^, sin 2 2<9 23 > 0.92, 
and sin 2 201% < 0.19, 5 ) with the exception of the LSND experiment, 17 ) whose re- 
sults were not confirmed by either the KARMEN experiment 18 ' or the MiniBooNE 
experiment. 19 ^ We do not take account of the ambiguities of these parameters in 
this paper for the sake of a clear presentation of our idea. The influence of these 
ambiguities will be discussed in 

We consider the search for the CP- 
violating phase 6 and the mass hierarchy 
sgnJm^ by investigating the — > v e 
appearance probability in long baseline 
neutrino oscillation experiments. We 
do not treat the disappearance chan- 
nel, which is effective in practical anal- 
yses to restrict parameter values such 
as the absolute value of dm^. Let us 
first illustrate with an example that the 
— ► v e appearance probability enables 
to search the value of 5 and the sign of 
8m\]_. In Fig. [JJ we display the — > u e 
appearance probabilities in a case of a 
baseline length of L = 700 km. There, 
the value of 5 and the sign of Sm^ are 
varied while other parameters are fixed 
to a set of example values: 
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The matter density p on the baseline is 
assumed to be constant and fixed to 



Neutrino Energy / [GeV] 

Fig. 1. The —> v a appearance probabilities 
and the trajectories of their first peaks for 
a baseline length of 700 km. The parame- 
ter values in Eq. (|2-1[) are adopted. Figure 
(a) is for the normal hierarchy and Fig.(b) 
is for the inverted. Each includes the prob- 
ability spectra for 8 = 0,7r/2,7r and 37r/2. 
Figure (c) overlays the trajectories of the 
top two, along with two oscillation spectra 
peaked at the intersection of these trajec- 
tories marked by the "x" symbol. 



/ o = 2.6g/cm d , (2- If) 

which is related to the electron num- 
ber density n e as n e = NxY e p with the 

Avogadro constant N\ and the proton-to-nucleon ratio Y e on the baseline. For con- 
venience, we use this set of values and assume Y e = 0.5 for the numerical calculations 
in this paper unless otherwise noted. Figure [JJ plots the energy spectra of the appear- 
ance probability (a) for dm^ = +2.5 x 10 -3 eV 2 > (normal hierarchy) and (b) for 
Sm^ = —2.5 x 10~ 3 eV 2 < (inverted hierarchy). Figure QJc) displays two spectra 
for the set of parameter values given in it. We can clearly see the dependence of 
the spectrum upon 5 and sgn 5m\i in (a) and (b), and thus we can search for these 
values through — > v e appearance experiments. 
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This dependence, however, does not guarantee that we can uniquely determine 
these values from a single experiment. Experimental results in some cases are con- 
sistent with both normal and inverted hierarchies, each with an appropriate choice 
of the value of S. For example, an experiment providing the oscillation spectrum 
of the dashed curve in Fig. QJc) would reject neither (sgn^m^, <5) = (+,0.637r) nor 
(sgn 8m\i , 5) = (— ,1.98-71") and leads to the hierarchy degeneracy, if the low-energy 
neutrinos below about 1 GeV are unobservable. The degeneracy introduces unfavor- 
able complication into the analysis and should be avoided. 

We now direct our attention to the peak of the appearance probability spectrum 
since it gives insight into the presence of hierarchy degeneracy. An appearance 
probability has a series of peaks (£ pea k,n, Ppeak,n) (n = 0, 1,2, •••) with E pcak ,o > 
E p eak,i > -E-peak,2 > ' ' ' • Each peak traces out a closed trajectory as we vary the 
value of 5 from to 2ir, keeping the other parameters fixed. The loops traced by 
the first peak (n = 0) are presented in Fig. [2(a) for dm^ > and (b) for Sm^ < 0. 
The peak positions for 5 = 0, tt/2, tt, and 37r/2 are indicated on the loops by 
open triangles, solid triangles, open diamonds, and solid diamonds, respectively. We 
observe that the peak position for the normal (inverted) hierarchy moves clockwise 
(counterclockwise) on the loop as the value of 5 increases. 

In Fig. [T](c), the peak loops for Sm^ ^ are copied from Figs. [2(a) and (b) with 
the symbol "x" marking one of their intersections. This intersection corresponds 
to the two sets of parameter values given in the figure. With these values, the 
shown oscillation probabilities are both peaked at the symbol as should be the case. 
We observe the similarity of the two for E > 1 GeV despite their difference for lower 
energies. Given the typical visible energy of neutrinos E > (0.5-1.0) GeV, we expect 
that their similarity will make it difficult to distinguish the two sets of parameter 
values by experiments and thus lead to degeneracy. 

We arrive here at an intuitive understanding of the presence of degeneracy. 
Assume that the visible energy range of an experiment covers a peak (-E p eak,n> -Ppeak.n) 
of the appearance probability. (In £j5]we discuss cases in which two or more peaks 
are visible, which occurs when a sufficiently long baseline, typically L > 1000 km, 
and a wide-band neutrino beam are used.) The values of E pe& ^ n and -P pea k,n depend 
on the oscillation parameters such as 5mf-, 9ij, and 5, which we collectively denote 
by {$i}- The observation made in the previous paragraph signifies that the two 
parameter sets and will be degenerate when the peak-matching condition 



is satisfied. This condition was first tested in Ref. 10) in the analysis of an example 
study, although not extensively. In particular, the hierarchy degeneracy arises when 
the values of Sm^ in and in have opposite signs. 

We confirm the validity of the intuitive discussion above by carrying out a quan- 
titative comparison of the two spectra shown in Fig. [2(c), following the analysis 
given in Sec. Ill A of Ref. 10). Our analysis assumes the neutrino beam flux shown 
in Fig. [21 which is the spectrum of the wide-band neutrino beam expected from the 
upgraded Alternate Gradient Synchrotron at Brookhaven National Laboratory in 



-Epeak,n = -Epeak,n({$i}) , 

({*}) =Ppeak,n({<}) 



(2-2a) 
(2-2b) 
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Fig. 2. The neutrino fluxes used in the exam- 
ple analysis. The flux of muon neutrinos 
(solid curve) is the spectrum of the wide- 
band neutrino beam expected from the up- 
graded Alternate Gradient Synchrotron at 
Brookhaven National Laboratory. 20 - ) The 
flux of electron neutrinos (dotted curve) is 
assumed to be 0.7% of that of muon neu- 
trinos with the same energy dependence. 
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Fig. 3. Allowed regions for 68.3% (gray curve) 
and 95% (black curve) confidence levels 
obtained from an example \ 2 analysis of 
the z/ M — > u c appearance events between 
0.7 GeV and 4.1 GeV. The large "x" sym- 
bol in the top figure indicates the true 
parameter values of our choice and the 
small one in the bottom indicates the peak- 
matching partner of the true values as ex- 
plained in the text. 



the U.S.A. 20 ) This flux is suitable for the present example where we continue to use 
L = 700 km, since its maximal energy E ~ (1~2) GeV coincides with the first-peak 
energy of the oscillation probability at this distance for l^m,^ | ~ (2-3) x 10~ 3 eV 2 . 
We also assume a water Cerenkov detector with a fiducial mass of 500 kt and a data 
acquisition time of 5 x 10 7 s. Other details of the setups are the same as in Ref. 10). 
The analysis is briefly outlined as follows. We fix the set of "true" parameter values 
as {drn^i, S) = (+2.5 x 10~ 3 eV 2 , 0.637r) taken from Fig. QJc), along with the values in 
Eq. (JED). The "test" values of 6 and Sm^ are varied over the parameter space while 
other parameters are fixed to their true values. We generate an energy spectrum for 
this test value and check its consistency with the true value via a x 2 goodness-of-fit 
analysis. The test values that pass this check constitute the allowed region. Figure 
[3] shows the allowed regions with 68.3% and 95% confidence levels obtained from the 
X 2 analysis over the range 0.7 GeV < E < 4.1 GeV. The top graph is for the normal 
hierarchy, and the bottom for the inverted. The large "x" symbol on the top graph 
indicates the true values of the parameters, and the allowed region extends around 
it due to statistical error, systematic error, and the background. Also in the bottom 
graph is the allowed region, which indicates the presence of a hierarchy degeneracy. 
The small "x" symbol located at [5mlx, 8) = (-2.5 x 10~ 3 eV 2 , -0.02vr) indicates 
the peak-matching partner of the true value of our choice shown in Fig. HJ^c) as 
— 0.027T = 1.987T mod 2tt. Note that it falls just within the allowed region with the 
wrong hierarchy, indicating the degeneracy of the two sets of parameter values. 
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§3. Degeneracy in light of the peak loops 



We showed in the previous section that the degeneracy follows from the two 
parameter sets which bring the oscillation peak to the same position. We employ 
the peak loops in this section to present a simple understanding of the presence of 
the degeneracy. 



3.1. Emergence of degeneracy: case studies 
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Fig. 4. Trajectories of the first oscillation 
peak for the baseline length of 700 km. 
The oscillation parameters in Eq. (|2-1[) are 
adopted. The solid curve is for the normal 
hierarchy and the dotted one for the in- 
verted. The "x" and "+" symbols denote 
the intersections of the trajectories. The 
points of A, A', B, B', C, and C' are re- 
ferred to in the text. 



We explain how we can read the 
presence and absence of the hierarchy 
degeneracy from the peak loops intro- 
duced in the previous section. In Fig. |U 
we reproduce the pair of peak loops 
presented in Fig. CD!c). The upper 
loop for 5m\i > and the lower loop 
for 8m%i < intersect at two points, 
each of which gives a pair of values of 
(sgn Srri^i, <5); in this example the two 
pairs are 



{(+,0.63tt),(- 1.98tt)} , 
{(+,0.127r),(-,1.34vr)} . 



(3-1) 



The true values of (sgn Sm^ , 5) pro- 
vided by Nature give an oscillation spec- 
trum whose first peak falls on some 
point on the two loops. We can cor- 
rectly determine these values by iden- 
tifying this point through experiments. 
Sufficient precision and accuracy are necessary to put the determination into prac- 
tice. Even well-controlled experiment, however, can fail to determine uniquely the 
parameter values if the true peak falls right upon an intersection of the loops, giving 
two possible sets of values as seen in Eq. ()3Tp . We demonstrate how degeneracies 
emerge in three typical cases where the true peak is at the points A, B, and C in 
Fig.H 

Case (A): The true hierarchy is normal and the true oscillation peak is at the 
point A in Fig. [H Owing to the vertical separation of the two loops, the sign of Sm^ 
is determined as positive once the appearance probability is determined precisely 
enough to distinguish the point A from the point A'. The value of 5 is restricted to 
a single allowed region around the value for A and the size of this region depends on 
the experimental precision and accuracy. 

Case (B): The true hierarchy is normal and the true oscillation peak is at the 
point B in Fig. |U The value of 5 is constrained to a single allowed region as in the 
previous case if the energy resolution of the experiment is high enough to distinguish 
the points B and B' on the parallel sides of the loop. If, on the other hand, the 
resolution is not sufficient, the values of S for B and for B' become indistinguishable 
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and the allowed region will extend around these two values. 

Case (C): The true hierarchy is normal, and the true oscillation peak is at the 
point C, which is the intersection of the two loops. There exists another value of 5 
that brings the oscillation peak to the same point C but with the inverted hierarchy, 
and we are thus led to a hierarchy degeneracy. Here we regard the points C and 
C to be on the upper and lower loops, respectively, in spite of their coincident 
positions. The allowed region consists of two separate parts which extend around 
the parameters for C and for C. 

3.2. Degeneracy in varying the baseline length 

Next we examine how a pair of loops 
move and become distorted as we vary 
the baseline length of the neutrino oscil- 
lation experiments. We present in Fig. [5] 
peak loops obtained as we vary the 
baseline length from 300 km to 1500 km, 
while fixing the other parameters to the 
example values given in Eq. (|2-ip . For 
a relatively short baseline (L ~ 300km 
in our example), the loop for the nor- 
mal hierarchy (the "normal loop") and 
for the inverted hierarchy (the "inverted 
loop") have similar shapes and overlap 
each other significantly. As the baseline 
becomes longer, the normal loop and the 
inverted loop move upward and down- 
ward, respectively, and they diverge due 
to the matter effect. The inverted loop 
is at the same time appreciably stretched in the Indirection and flattened in the 
P-direction. The two loops are seen to become disjoint at a certain baseline length. 
This critical length L C rit for the present parameter set is L cr it — (1100- 1300) km as 
can be read in Fig. 

From the above observations, we reach following outlook upon the determination 
of the hierarchy and the CP-violating phase. 

The hierarchy is difficult to determine by experiments with a short baseline 
length where the pair of loops overlaps considerably. We expect, however, that 
the hierarchy can be determined when the true peak is located at the top part of 
the upper loop or at the bottom part of the lower loop. The determination of the 
hierarchy becomes easier as the baseline becomes longer and as the increasing effect 
of the matter separates the pair of loops. When the baseline is longer than the critical 
length, the condition given in Eq. (I2-2D is never satisfied with opposite hierarchies 
and the hierarchy can thus be uniquely determined regardless of the value of S. 

The search for the value of 5 is entangled with that for the hierarchy. The 
hierarchy degeneracy is one obstacle to the determination of the value of 5 especially 
when the baseline is short. Another obstacle is the experimental limitation on the 




Peak Energy / [GeV] 

Fig. 5. Trajectories of the first oscillation 
peak for baseline lengths between 300 km 
and 1500 km. The oscillation parameters 
in Eq. (|2T|) are adopted. The solid curve 
is for the normal hierarchy and the dotted 
one for the inverted. 
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energy resolution and on the precision and the accuracy of the oscillation probability. 
The loop is narrow in the indirection when the baseline is short, and the resolution 
of the energy must be sufficiently high to distinguish the parallel sides of the loops; 
otherwise an additional degeneracy will be introduced. In contrast, the inverted loop 
is flattened in the P-direction when the baseline is long. We thus require precise and 
accurate measurements of the oscillation probability to avoid introducing an extra 
degeneracy. 



§4. Analytic study of the peak loops 



We have illustrated the determination of parameters with the aid of peak loops 
and found their position, size, and shape informative. In this section, we exploit 
analytic expressions of the oscillation probability to formulate them in terms of 
the mass parameters, the mixing parameters, and the baseline length. We thereby 
analyze how a pair of peak loops for Sm^ ^ are distorted and separated as the 
baseline length increases. 



4.1. Formulation of the peak loops via the oscillation probability formula 

We derive formulae for the oscillation probability with a perturbative expansion 
of the 5-matrix in terms of A 2 \ = 5m 21 L/2E and A m = \^2Gpn c L, where Gp = 
1.166 x 1CP 5 GeV~ 2 is the Fermi constant. The derivation was first worked out to 
first order to analyze the effect of CP violation separately from the matter effect. 21 ) 



In the present paper, we calculate the Vn 



appearance probability to second 



order with the following consideration. Note that a peak loop collapses to a point 
when 5m 21 vanishes since CP violation is a three- generation effect. The size of the 
peak loop is thus of first order, and its distortion depending on the baseline length, 
in which we are interested, is of second order or higher. The second-order calculation 
outlined in Appendix [A] results in lengthy expressions and we apply an additional 
simplification that exploits the smallness of #13. We drop O (sin 2 #13) terms in the 
coefficients of A 2l and A m A 2 i, as well as 0(sin 3 #i3) terms in that of A^. Here we 
take account of the relation A m > A21, which holds for the cases in which we are 
interested (see Appendix [B]). We then obtain 



Pk -> u e ,E) = 4l(Asm 2 + B) , 



(4-la) 
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and 



B 



^ll J 2 • 2, 

— ^sm 5. 



(4-ld) 



Here A 



6mlL/2E, I = c?,s? " 2 



-13*13*23) 



S23, and j = C13S13C23S23C12S12, where s. 



sinl 



and Cjj = cos 0y . The approximation we employed in deriving Eq. (|4-ip is suitable 
especially for a short baseline, typically for L < O(10 3 km); see Appendix [Bl The 
oscillation probability for anti-neutrinos P{pn — * v e ,E) is obtained by changing the 
signs of 5 and A m . 

We obtain the peak energy, which gives the local maxima of the oscillation 
probability of Eq. (|4-1|) , as 
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where n = 0, 1, 2, • • • is the peak index, 77 = (2n + l)ir, R = Sm?,i/\dm^i\, and the 
top of the double sign is for Sm^i > and the bottom for Sm^ < 0. The oscillation 
probability at the peak energy is given by 
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(4-3) 



A peak loop is obtained when we keep track of (E pea ^ >n , P pe ak,n) as we vary 5 with 
other parameters fixed. We present in Fig. [U]the peak loops for the first peak (n = 0) 
obtained from Eqs. (|4-2p and (|4 3j) and compare them with the numerical results to 
confirm good agreement. 

The basic properties of the peak loop are its position, size, and shape, and we 
quantify them using its central position, width, and height, respectively, as follows. 
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Fig. 6. Trajectories of the first oscillation 
peak for baseline lengths between 300 km 
and 1500 km obtained from the numerical 
calculation (the solid curve for the normal 
hierarchy and the dotted for the inverted) 
and analytic approximations (the dashed 
curve for the normal hierarchy and the 
dash-dotted for the inverted). The oscil- 
lation parameters in Eq. (|2-1[1 are adopted. 




Fig. 7. Notations to designate the central po- 
sition, the width, and the height of the tra- 
jectories of the (n + l)-th oscillation peak. 



The terms independent of 5 in the braces of Eqs. (|4-2|) and (|4 3j) give the average 
values of -E p0 ak,n and -P pca k,n over a cycle of 5 and can be regarded as the center of 
a peak loop. The central position up to first order thus reads 



l 77 



4/ 



l±A m -{l-2s( 3 ) 



12 



(4-4) 



where the double sign is the same as in Eqs. (|4-2p and (|4-3p ; see Fig. [JJ The terms 
that depend on 5 in Eqs. (|4-2p and (|4-3|) account for the size and the shape of a 



peak loop. The width of a loop AE^ and its height AP^ can be estimated by 
taking the difference of the maximum and minimum values of -Ep ea k,n an d 7p e ak,n as 
functions of 5 (see Fig. [7J . They are calculated up to second order as 



AE^ =\5m 2 31 \LR^\l 1 + 



1 

77 



1 T A 



2 1 



77 2 I 
32/ n 4 



77 1 + 4/77 2 1 + 4/77 2 



Z\PW = 87?77j 



l±Z\ m -(l 



(4-5a) 



(4-5b) 



We remark on the magnitude of the correction terms appearing in Eqs. (|4-4p 
and (|4-5p . The correction terms include the matter effect proportional to A m and 
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the three-generation effect proportional to R. Of these two, the three-generation 
effect is even smaller than the matter effect for a typical long baseline experiments. 
Evaluation using the example values given in Eq. (|2-ip is sufficient to verify this with 
an order-of-magnitude estimation, which yields 



l 



0.2 



L 



1000 [km] 



for n = , R(l - 2s\ 2 ) ~ Rs\ 2 ~ 0.01 . 



in 



2s j 3 ) T 



4/ 



1 + Z\ m -(1-2 S 2 3 ) 



(4-6) 



The expressions for anti-neutrinos are obtained from Eqs. (|4-4p and (|4-5p by 
simply flipping the signs of A m : 



(4-7) 



AE^^mi^LR^l + jpl 
2 1 - 32/ II 4 



1± A 



zAP^ = SRIIj 



1 + 



77 1 + 4/77 2 ^ 1 + 4/77 2 
A. — fl-— 1 ±Ps1 



77 V 77V 



(4-8a) 



(4-8b) 



Equations (|4-7p and (|4-8p have an interesting relation to (|4-4p and (|4-5p . The former 
pair for anti-neutrinos is equivalent to the latter for neutrinos with flipped hierarchy 



up to the correction terms proportional to Rs\ 2 or R{\ 
shown in Eq. (|4-6p . 



2s 



127' 



which are small as 



4.2. Analysis of the peak loops 

We elevate the observation of peak loops for neutrinos shown by Fig. [5] to a sys- 
tematic analysis by applying the formulae we developed in the previous subsection. 
We also extend our consideration to the anti-neutrino case. 

Equation (|4-4p reveals the order behind the motion of the loops for neutrinos 
in the E-P plane as the baseline length becomes longer. The leading dependence 
of on the baseline length comes from the prefactor \5m 31 \L/2E and drives the 
loops rightward. Its subleading dependence due to the matter-effect term, which is 
proportional to the baseline length through A m and to sgn 6m 31 , provides a correc- 
tion which pulls the normal loop back leftward and gives the inverted loop another 
push rightward. The other subleading term Rs\ 2 of E^ 1 gives an extra correction 
whose sign depends on the hierarchy and whose magnitude depends on the mass 
parameters and the mixing parameters but not on the baseline length. The depen- 
dence of on the baseline length comes from the matter-effect term, which raises 
the normal loop and lowers the inverted loop. As a whole, a pair of normal and 
inverted loops are driven rightward while they split vertically, and their alignment 
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tilts counterclockwise. These features of the motion of the loops are clearly observed 
in Fig. [5j 

We next analyze the distortion of 
the loop for neutrinos in terms of its 
width and height given in Eq. (|4-5p . The 

expression for AE^ has an overall fac- 
tor which is proportional to the base- 
line length and accounts for the widen- 
ing of the loops. One of two sublead- 
ing corrections of AE^ is the matter- 
effect term, which adds a dependence on 
the baseline length. The other term also 
gives a correction, but does not give any 
additional dependence on the baseline 
length. Both of these terms depend on 
the hierarchy, and decelerate the widen- 
ing of the normal loop and accelerate 
that of the inverted loop. The depen- 
dence of AP^ on the baseline length 
is due to the subleading matter-effect 
term, which vertically stretches the normal loop and compresses the inverted loop, 
and the other correction term supplements its effect. All these features of the loops 
are evident in Fig. [5j 

Now that the motion and the distortion of the loops for the neutrinos have 
been analyzed, those for anti-neutrinos are easy to derive owing to the discussion 
below Eq. (|4-8|) in the previous subsection. A loop for anti-neutrinos moves with 
distortion just as that for neutrinos with the opposite hierarchy, except for the small 
contribution of the three-generation effect. 

The critical baseline length mentioned in the previous section also can be ana- 
lyzed with our formulation. The critical length is defined as the maximum length 
giving the intersection of the pairing loops with opposite hierarchies. It is calculated 
in the first-order approximation as 



Fig. 8. The critical length as a function of 
sin 2 2#i3 in the first-order approximation, 
Eq. (|4-9|) . The solid curve is for neutri- 
nos and the dotted curve for anti-neutrinos. 
We adopt the values of the oscillation pa- 
rameters given in Eq. (|2TI) except for that 
of sin 2 2fli3 in Eq. (l2Te|) . 



L, 



1 



n 



crit 



V2G 



F n e 



12 64 

i72 + 7J4 



R 



C23C12S12 



1 



S13S23 1 



?s 2 



12 64 

lP + W 



4 / S13S23S12 



Tin 



C23C12 



T 1 



_8_\ £13£23£12 

n 2 ) C23C12 



(4-9) 



where the double sign is — for neutrinos and + for anti-neutrinos. The prefactor ap- 
pearing on the right-hand side of this expression gives \/\/2G-pn c = 5.17 x 10 3 [km] x 
(p/[g cm -3 ]) -1 = 2.0 x 10 3 [km] • (p/[2.6g cm -3 ]) -1 . The critical length is inversely 
proportional to sin #13 owing to the factor C23C12S12/S13S23, with the small correction 
of 0(^13). It is particularly important that this dependence can make the critical 
length very long, since the lower bound on the value of sin #13 is unknown at present. 
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Figure [8] shows the dependence of the critical length on sin 2 2#i3 for the first peak, 
where the other parameters are fixed to the values given in Eq. (|2Tp . The critical 
length for anti-neutrinos is slightly longer than that for neutrinos due to the small 
correction proportional to si3S23.si2/c23Ci2. From this graph, we find -L cr it = 1150 km 
for neutrinos and L C rit = 1200 km for anti-neutrinos at sin 2 2#i3 = 0.06. This graph 
is qualitatively consistent with Fig. 5 in Ref. 16), which is obtained from their bi- 
channel analysis with the energy fixed around a peak at E = \8m^y\L/2ir. 

It is remarkable that the following straightforward analysis of the loops repro- 
duces the main part of Eq. (|4-9fl . which is tedious to derive, and clarifies the origin 
of the dependence of the critical length upon the mass parameters and the mixing 
parameters. We compare the separation between the loops with their size to obtain 
a condition for their disentanglement. The separation of the two loops is estimated 
using the difference between the respective values of E and P , while the size 
of each loop is measured by its width and height. We consider the ratio of the sep- 
aration to the size to compare the two, and define a normalized separation vector 
as 



[AE ( n +) + AE { n ] ] /2 ' [APi +) + APt ] ] /2 

-^(1-2,?,/ 1 l ~ ^ 



R j V n y/i + 4/Ti 2 ' n 2 



(4-10) 



where we have omitted small terms such as i?s 2 2 and R(l — s 2 2 ) on the right-hand 
side for simplicity. The magnitude of this vector is given by 

which is proportional to L and parametrizes the configuration of the paired loops, 
i.e. intersected, tangential, or disjoint. The normalization of Eq. (|4T0p is chosen so 
that the separation at the critical length S CI it = S(L CT n) becomes unity when the 
loops are two similar ellipses whose major axes are parallel. We write the critical 
length from Eq. (|4- as 

(2n + l)7T C23C12S12 1 1 q , 
-^crit — — 7=7; -K ~ '-'crit- \^' ll > 

V2G F n e s 23 si3 1 - 2sf 3 

Equation (|4- 12|) approximately reproduces the dependence on the parameters ap- 
pearing on the right-hand side of Eq. (|4-9p . assuming that 5 cr it does not give any 
extra dependence on the mass parameters and the mixing parameters. 

§5. Toward a resolution of the hierarchy degeneracy 



In this section, we consider how to avoid the hierarchy degeneracy in determining 
the value of 5 with long baseline experiments. 
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Fig. 9 



A promising approach for determin- 
ing the value of 5 and the hierarchy is (1) 
to carry out an experiment with baseline 
length longer than the critical length so 
that the hierarchy is identified regard- 
less of the value of 5. This approach 
can be made feasible by pinning down 
the position of only one peak. It is 
nonetheless not free of hurdles: the re- 
quired baseline length is typically about 
1000 km and may be even longer, de- 
pending on the value of #13. Perform- 
ing experiments with such a long base- 
line is challenging due to, for instance, 
the small flux of the neutrino beam and 
a possible large ambiguity of matter ef- 
fects. 

Experiments with shorter baselines 
are more feasible, but leave the possi- 
bility of degeneracy. To overcome this 
problem, we can simultaneously make use of two or more peaks. We examine the 
effectiveness of the following three approaches of this kind: (2) observing — > v G 



Trajectories of the first peak of (a) the 
— » v c probability and (b) the —* v c 
probability for a baseline length of 700 km. 
The solid and dotted curves are for the nor- 
mal and inverted hierarchies, respectively. 
The "x" and "+" symbols in (a) are on 
the intersections of the peak loops, and 
the corresponding values of the parameters 
are also shown. The symbols in (b) mark 
the points corresponding to these parame- 
ter values. 



appearance events in addition to 



events; (3) carrying out two or more ex- 



periments with different baseline lengths; and (4) carrying out an experiment which 
has two or more oscillation peaks within its range of visible neutrino energy. 

The second approach employs both neutrinos and anti- neutrinos. We analyze 
the merit of this approach with Fig. [9l which presents two pairs of peak loops, one 
(a) for neutrinos and the other (b) for anti-neutrinos, with the values of Eq. (|2-ip 
and L = 700 km. The loops in Figs. [9^a) and (b) resemble each other with opposite 
assignments of the hierarchy, confirming our discussion in §4.21 based on analytic 
formulae. Marked by the "x" and "+" symbols in Fig.G^a) are the two intersections 
of the normal and inverted loops. Each intersection is associated with two sets of 
(sgn 6171^,5) as we gave examples in Eq. (|3-ip and corresponds to the presence of the 
hierarchy degeneracy. Turning to Fig.[9](b), the four parameter sets considered there 
correspond to four distinct points on the peak loops for anti-neutrinos. Hence the 
combined analysis of — > v e events and 9^ — ► u e events will be able to resolve the 
hierarchy degeneracy, provided that the two "x" symbols or the two "+" symbols 
are experimentally distinguishable. In our present example, the two "+" symbols 
for anti-neutrinos are close to each other in the E-P plane and are thus difficult to 
distinguish. If indistinguishable, the hierarchy degeneracy will persist even with the 
aid of anti-neutrino events. 

The third approach makes use of two or more different baseline lengths. 15 ^ We 
reproduce in Fig. [10] the series of peak loops shown in Fig. On top of them, 
we mark the points corresponding to the parameter values for the intersections of 
the loops in the case L = 700 km as we did in Fig. [9l Each pair of loops for a 
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(sgn(8my, 8c) = 
X ((+,0.63l8,(— , 1.981)] 
+ ((+, 0.12K), (—, 1.34i)l 



■ Sm?,>0 

■ Smi, < 

5 C p = 
a 
A 0.571 



1300 1500 km 



12 3 4 

Peak Energy / [GeV] 

Fig. 10. The positions of the first peak 
for baseline lengths between 300 km and 
1500km. The "x" and "+" symbols in- 
dicate the positions for the four sets of 
parameter values which correspond to the 
two intersections of the peak loops at L — 
700 km, and they are plotted on top of the 
series of peak loops reproduced from Fig. [5] 
The parameter values given in Eq. (|2T[1 are 
adopted. The solid (dotted) curve is for the 
normal (inverted) hierarchy. 
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Fig. 11. Trajectories of the first and second 
oscillation peaks for L — 1100 km. The 
parameter values given in Eq. (|2T[) are 
adopted. The solid (dotted) loop is for the 
normal (inverted) hierarchy. The "x" and 
"+" symbols correspond to the sets of pa- 
rameter values for which the first peak is 
at the intersection of its trajectories. The 
points for these values are also on the sec- 
ond peak loops and are marked by the cor- 
responding symbols. The oscillation spec- 
tra for the parameter values are overlaid. 



fixed baseline length has two "x" symbols and two "+" symbols, which reduce, by 
definition, to a single "x" symbol and a single "+" symbol at L = 700 km. The 
separation between the two "x" symbols as well as between the two "+" symbols 
grows as the baseline length deviates from 700km. The degeneracy of the two "x" 
symbols or the two "+" symbols can thus be resolved with an additional experiment 
with a different baseline length, where the large difference between the two baseline 
lengths is preferable. 

The fourth approach exploits other peaks of the appearance probability in ad- 
dition to the first one for n = 0. 20 ^ This approach requires a wide-band neutrino 
beam in order to accumulate a sufficient number of events at the multiple peaks. An 
analysis covering peaks for different values of n is effective to resolve the degeneracy 
since each peak has its own properties due to the factor 77 = (2n + 1)tt in Eqs. (|4-4p 
and (|4-5p . The peak for n > 0, however, is not easy to observe due to its small energy 
which is suppressed by the factor of l/(2n + 1). On the other hand, elongation of the 
baseline makes the peak energies higher, bringing both the first and second peaks 
within the range of visible energy. We show in Fig. [TT] the two sets of loops traced 
by the first and second peaks for L = 1100 km. The second-peak loops are stretched 
in the P-direction and have a significant overlap, showing their noticeable difference 
from the first-peak loops. The "x" and "+" symbols correspond to the values at 
the two intersections of the loops for the first peak. On the second peak loops, the 
paired "x" symbols as well as the paired "+" symbol are separated from each other 
and are experimentally distinguishable. The hierarchy degeneracy in the first peak 
will be thereby removed by observing the second peak. 
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§6. Conclusion and discussion 

We have studied the search for the leptonic CP- violating phase 5 and the neu- 
trino mass hierarchy sgndm^ through the observation of the — > v c oscillation 
with long baseline experiments. The energy spectrum of the v„ — > u e appearance 
probability is engraved with the values of these parameters, although not uniquely 
The search for these parameters through this spectrum can lead to degeneracy when 
the two spectra corresponding to the two different parameter values are indistin- 
guishable. We implemented the presence of degeneracy using the following criterion: 
the oscillation probabilities for the two parameter values are peaked at the same en- 
ergy and have the same peak probability. In light of this peak-matching condition, 
we systematically examined the presence and absence of degeneracy, especially the 
hierarchy degeneracy with regard to the value of 5 and the sign of Sm^i ■ We have also 
demonstrated the prospects for resolving the degeneracy using a single experiment 
and combinations of experiments. 

We introduced a looped trajectory traced by the position of the oscillation peak 
in the E-P plane while varying the value of 5 over the interval [0, 2tt]. We have 
shown that the loop plays a key role in gaining an intuitive understanding of the 
parameter degeneracy. We drew pairs of loops for both hierarchies and investigated 
their properties and behavior to obtain an overview of the emergence of the hierarchy 
degeneracy and methods to avoid it. We observed that paired loops with different 
hierarchies are completely separated when the baseline is longer than a critical length, 
which is proportional to 1/ sin #13 and is typically about 1000 km or longer. The 
hierarchy degeneracy can be avoided in experiments with a baseline length longer 
than the critical length. 

Employing the peak loops, we obtained a general understanding of experimental 
capabilities to determine the value of 5 and the mass hierarchy. We described four 
ideas for resolving degeneracy with long baseline experiments: using a baseline longer 
than the critical length, employing anti-neutrinos as well as neutrinos, using two or 
more different baseline lengths, and carrying out experiments covering two or more 
oscillation peaks. We discussed the merits and demerits of each idea in terms of the 
peak loops. We pointed out that the first idea of using a sufficiently long baseline is 
simple but is not free of hurdles. We also found for the second case that employing 
anti-neutrinos does not always resolve the hierarchy degeneracy. 

The peak loops we introduced to analyze the parameter degeneracy is as versa- 
tile as the trajectories in the bi-channel plots presented in Ref. 12). The two have 
differences at the same time. One difference is that our plot employs only a single 
channel and thus can be used to explore the capabilities of single-channel experi- 
ments. Another difference is that we make essential use of the oscillation peak as 
a representative of the oscillation spectrum over a finite energy range, while the bi- 
channel plot is drawn for an arbitrary fixed energy or for integrated values over an 
energy range at one's convenience. The significance of the peak position consists in 
its implication to the values of the parameters for which we search. Supported by the 
peak-matching condition, our method facilitates a direct prediction of the presence 
of degeneracy. 
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We kept the values of the mass parameters and the mixing parameters fixed 
except for 5 and sgn Sm^ , assuming that their values will be determined before the 
experiments in question are carried out. It is possible, however, that they will not 
be determined with sufficient precision by that time and that their ambiguities will 
add another complication to the searches. 9 )> 10 ) The peak loops presented in this 
paper also offer an intuitive understanding of the effects of these ambiguities. The 
position of the oscillation peak moves around in the E-P plane when the value of a 
parameter is ambiguous and varies within its allowed region. The peak position will 
thus be smeared and consequently the loop of the peak will appear broadened and 
blurred. The value of S will be accordingly obscured and the mass hierarchy will 
be misidentified within the extended region around the intersection of the loops. In 
addition, the broadening of a narrow loop makes its parallel sides indistinguishable. 
We hence expect that ambiguities of the mass parameters and the mixing param- 
eters will worsen the degeneracy of the target parameters and hinder attempts to 
determine them. We need extra effort to keep the ambiguities under control, such as 
combining with the — ► disappearance channel or employing reactor neutrino 
experiments. Detailed analyses of the ambiguities in the framework of the peak loops 
are left for future works. 

The present analysis based on the peak loops for the oscillation probability is 
independent of experimental specifications and gives an organized understanding of 
the presence and absence of degeneracy. A more practical evaluation of quantitative 
capabilities can be carried out for specific experiments by applying the expected 
number of events to our analysis in place of the oscillation probability. For that 
purpose, we need to calculate the expected number of events by using the neutrino 
beam flux of the experiment, the neutrino cross sections, the detector design in- 
cluding a knowledge of systematic errors and background, and other aspects of the 
experimental setup. Drawing the loops traced by the peaks of the spectrum of the 
event number, we can study the possibility of the presence of degeneracy, quantita- 
tively estimate the allowed region in the parameter space, and select an appropriate 
baseline length to determine the parameter values. 
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Appendix A 

Oscillation Probabilities to Second Order 



We calculate the S'-matrix of neutrino oscillation to second order in A21 = 
bmWLjlE and A m = aL/2E, and also derive a formula for the oscillation probability 
to this order. 

The evolution equation of neutrinos v = (v c , u^, v T ) T is given by 

i^=H{x)u{x), (AT) 
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where 



H{x) = ± (Vdiag(0, 6m 2 21 ,5m 2 31 )tf + dmg(a(x), 0, 0)) . (A-2) 

Here E is the neutrino energy, U is the unitary mixing matrix of the neutrinos 
(MNS matrix), and a(x) = 2\/2G-pn e {x)E represents the matter effect, where Gp 
is the Fermi constant and n e (x) is the number density of electrons on the baseline. 
The probability for a neutrino v a to change into up ({a, {3} C {e, /U,r}) is given by 



P(v a -> up) = \S{x)p a \ , 



where 



S(x) = T exp 



\J dsH(s)\ 



(A3) 



(A-4) 



We assume in this paper that the spatial variation of the electron density n e is 
negligible so that Eq. (|A-4p simplifies to 



-iHx 



S(x) = e 

We split H as H = Hq + Hi , where we have 

2^ C/diag(0 '°' 5m ' l)C/t ' 

Hl = 1e { Ud[ag (°> 5m n ' °) ui + diag(a ' °' 0) 

and treat Hi as a perturbation to obtain 



(A5) 

(A6) 
(A7) 



S(x) 



~ iHoX T exp (-i £ ds e iHoS Hie~ iHoS ^j . (A-? 



Expansion of Eq. (|A-8|) leads to a systematic approximation of S-matrix. 21 - 1 An 
explicit calculation of S(x) up to the second order in Hi gives 



where 



S\(x)p c 



S{x)j3 a = S Q (x)j3 a + S 1 (x)p a + 3 2 {x)p a , 

S (x)p a = S a p - Up 3 U* 3 {l - e~ iA ^) , 
iA 21 Up 2 U* 2 

SpeSae - (S ae + Sp e )Up 3 U* 3 + | U c3 1 2 Up 3 K 3 ( 1 + e"^ 31 ) 



A 

A 



(A9) 
(A10) 

(All) 



.31 



■(<5 Qe + <5/3e - 2|^ e3 | 2 )^3^3(l " e 31 ) 
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and 



S 2 {x)p a 

~ An^21 



0^21^/32^2 



-(5 ae + 6p e )Up 2 U ( 



<i2 



(U e2 U^Up 3 U* 2 + U^Up^) [\+ ' 



1 - e -^ 31 



An I g^oe ~ (o~ac + 6p e )Up 3 U* 3 ^-+^ 



1 



2 An 

% 1 - e~ izi31 



A 2 

^31 



All 



\Ues\ S - Sp e 5ae[ 77 + 



1 j 1 - e 



,2 Ai 
+ (1 + 5 Qe + Sp e )Up 3 U^ 



A 2 

^31 

1 i(2 + e~ iZi31 ) 3(l-e" izi31 

2 + a7 



A 2 

^31 



|C/e3| 4 ^3^3 



l_ e -^3i 3i(l + e- iZi31 ) 6(1 -e"^ 31 ) 



+ 



A 



31 



a 2 

^31 



(A-12) 



with A i:j = 5m 2 -L/2E and A m = aL/2E. Substituting Eqs. dAOO]) . (lA4T~|) and 
fA12l) into Eq. (IA-3jl . we obtain the expression for the oscillation probability up to 
second order as 



P(v a -* Up) = p (v a ~> Vp) + pi(u a -> l/g) + P2(^a ~> ^) , 



(A- 13a) 



where 



p (^ Q -> 1^3) = #a/3 - 4|LT /33 | 2 ((5 a/ 3 - |C/ a 3| 2 ) 



2 V. 



(A-13b) 



Pl(y a -» ^) = 4Z\ 2 i [Re {U*p 3 Up 2 U a3 U* a2 )sc - Im (U^ 3 Up 2 U a3 U* 2 

+ 4Z\ m |^e3| 2 (<5ae5 ( g e ~ ^oe 

- <M^ 3 | 2 + 2|f/ /33 | 2 |^ 3 | 2 ) f sc - -^-a 5 

V ^31 



* \„2 



(A-13c) 



and 



p 2 (>a -> ^) = A\\Vlk + Al^ 



(A-13d) 



with 



P2A = -(<5 a /3 - |^2| 2 )|t//3 2 | 2 + 2 Re (U*p 3 Up 2 U a3 U* 2 
+ 21m (U*p 3 Up 2 U a3 U* 2 )sc, 



(A-13e) 
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P2B = -2(5 ac + 6p e - 2\U e3 \ 2 ) Re {U* p3 U p2 U a3 U* a2 ) ( 1 - s 



A 31 



sc 



2S aP BB(U^U p2 UM)[l 



A 2 ' 

^31 



- 2 



- 2 



|^ 3 | 2 Re (U: 3 U e2 U a3 U* 2 ) + \U a3 \ 2 Re (U^U 02 U e3 U* e2/ 
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{5 ae + Spe - 2\U c3 \ 2 ) Im (U*p 3 U p2 U a3 U* 2 ) 



sc 



A 



31 



and 
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|f/e3| 2 (^e<5/3e " 5 ae \U(] 3 \ 2 - 5p e \U a3 \ 2 ) ( 1 - 2 S 2 

4 12 



Z\ 
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31 



12 

^31 



<5 ay g|?7e3| 1^/33 1 



1 + 



A 



-sc 



31 



Z\ 2 ' 
^31 



+ 2|J7 e3 | 2 |C/ /33 | 2 |^3| 2 M 



i. 
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-sc 



31 



12 

^31 



+ 2^|f/ e3 | 4 |^3| 2 ( S 2 + ^--c- Jj-s 2 ) 

+ |C/e3| 4 (<5ae^e " " Sp e \U a3 \ 2 ) U - 6s' 
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^31 
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Here we have used the abbreviations s = sin(Z\ 3 i/2) and c = cos(Z\ 3 i/2) 
The lengthy expression given in Eq. (|A-13d|) has the form 

T3( \ n ■ 2 Al n . Al Al . ^ 

P(i/ a -> i^) =Ci sin — - + C 2 sin — - cos — - + C 3 , 



(A-13f) 



(A-13g) 



(A-14) 



where the coefficients C±, C 2 , and C 3 are written in terms of the elements of the 
mixing matrix, the mass parameters, and A m . The coefficient C\ is O(l) while C 2 is 
first order or higher, so that Eq. (|A-14p reduces to a two-generation formula in the 
zeroth-order approximation. We factor out sin(Ai/2) from the sum of sinusoidal 
functions and transform the sum of the sine and cosine terms into a single sine term 
as 



A 



C'( + C'2 sin — ^- sin 



A 



31 



+ CA+ c 3 



C4 = arctan 



9i 

Ci : 



(A- 15a) 
(A-15b) 
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where we take < C4 < it for C2 > and it < C4 < 2ir for C 2 < 0. We rewrite 
Eq. (|A-15aP as 

P(v a - up) = y/(% + Ci sin 2 + 9±\ - + CI sin 2 9± + C 3 , (A-16) 

which reproduces the form of Eq. (|4-lap with 

A = i^T^,e=^ + ^,5 = I(c 3 -^T^sin^). (A-17) 

The expressions of Eqs. (f4~Ibl) . (jUTcJ), and dlTTd]) are obtained from Eq. fATTl) by 
setting a = fi and /? = e in Ci, C 2 , C3, and C4 and taking the terms up to the 
appropriate order. 

Appendix B 

Validity of the Approximation 



We made use of approximate formulae for the oscillation probability in £JH We 
elucidate the requirements for the these formulae to be valid. 
The approximation we developed in Appendix lAl requires 



5m\i 3> a,dm,2i , (B-l) 

An = — = 5.0 x ict 1 ! — -= — h( — — 5t)<i, (b-2) 

2E V 1000 t km ] / V 2 -6[gcm- 3 ] ' 



and 



A 21 (E) 



5m 21 L 
2E 

2,1 X 10_1 ^ ^-lO-sW 2 ]) (lOOO [km]) ( [GeV] ) <<; L 



(B-3) 



Since neutrino energy we consider in this paper is E Speaks ~ |^?Ti| 1 |-L/[2(2n + 
l)-7r], the conditions given in Eq. ([B]) are satisfied at least marginally for L ~ 
(300-1500)km. 

We additionally simplified the expression by dropping higher-order terms of 
sin #13 taking account of its smallness as mentioned in £J5J Taking account of the 
relation A 2 \ < A m , we dropped terms of 0(sin 2 #43) in the coefficients of A^i and of 
A 2 \A m , as well as terms of 0(sin 3 #43) in the coefficients of A\^. This simplification 
is valid when the omitted terms are much smaller than the leading term which is of 
O (sin 2 #13), or explicitly when the following conditions are satisfied: 

5ml 



l 21 



4ji(W^ ^ gag = L1 x 1Q -2 . (2n + 1)2 V8-2 • 10-5 [eV 2 ] ; ^ 



sin 2 6*43 / I <5m 2 



311 



2.5 • IO- 3 [eV 2 ] 
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4a(WAnsin'g 13 = 5 2 x 1Q _ a . (2n + 1) U2.10-5[eV*] 



5m 2 



21 



sin 2 6*i3 / | <5m 2 



31 



,2.5 -10- 3 [eV 2 ] 



(B-B) 



1000 [km]/ V2.6 [gem 



-31 



< 1 



and 

= 2.5 X 10- ■ sine 13 ( - X( / S « 1- P") 
sin 2 0i3 J Vl000[km]y V2-6 [g cm^ 3 ] J 

These inequalities are satisfied for allowed values of sin 2 2#i3 < 0.19. 5 ) 
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